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We present here a quantum mechanical framework for defining the statistics of measurements of 
J dtA(t), A(t) being a quantum mechanical variable. This is a generalization of the so-called full 
counting statistics proposed earlier for DC electric currents. 

We develop an influence functional formalism that allows us to study the quantum system along 
with the measuring device thus fully accounting for the action of the detector on the system to be 
measured. We define the full counting statistics of an arbitrary variable by means of an evolution 
operator that relates initial and final density matrices of the measuring device. 

In this way we are able to resolve inconsistencies that occur in earlier definitions. We suggest two 
i-^h , schemes whereby the so defined full statistics can be observed experimentally. 
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I. INTRODUCTION 



The measurement paradigm in quantum mechanics assumes at least implicitly that the measurement is instant .El 
This is to be contrasted with a realistic measurement of, say, electric current where the result of measurement is 
averaged over a sufficiently long time interval. If one intends to measure a variable A, the individual measurement 
gives L A(t)dt/T . The reason for this is obvious: any measurement has to be accurate. The integration over time 
reduces possible instant fluctuations of A(t) resulting in a more accurate outcome of an individual measurement of 
this sort. The dispersion of the probability distribution of the outcomes is supposed to vanish in the limit of T — > oo. 

•y-, , This paper focuses on the problems related to the determination of this probability distribution, the full statistics of 

£T} ' the measurement results. p n 

t-H ! Several years ago Levitov and Lesovikoa have made a significant step in the understanding of this fundamental 

issue. They have introduced a concept of full counting statistics (FCS) of electric current and have found this statistics 
for the generic case of a one-mode mesoscopic conductor. The word "counting" reflects the discreteness of the electric 
charge. If electrons were classical particles, one could just count electrons traversing a conductor. The full counting 
statistics could be readily defined in terms of a probability to have N electrons transferred through the conductor 
during a time interval T, Pr(N). With this distribution function one calculates the average current (N)/T, current 
noise ({N) 2 > —(N) 2 )/T and all higher cumulants of the current. A non-trivial value of interest is the probability 
to have big deviations from the average value. This can be measured with a threshold detector. The probability 

'■^J ■ distribution Pt(N) would be the goal of a quantummechanical calculation. 

The operator of electric current through a conductor, J,Js well-defined in the Fock space spanned by the scattering 
states of electrons. The initial idea of Lesovik and LevitovB was to define an operator of transferred charge by means 
of a seemingly obvious relation 



Qtr = / dtl(t). (1) 

Jo 

Having this operator in hands, one applies the general paradigm of quantum measurements: The probability to have 
a certain charge q transferred equals the square of the projection of the wave function of the system on the eigenstate 
of Q tr with eigenvalue q. Lesovik and Levitov were able to perform a challenging calculation of these projections. 
However, they were hardly satisfied with the results. For instance, the transferred charge was not quantized in units 
of the elementary charge. _ 

This is why in the subsequent paperlj the same authors have proposed another method of evaluating Pt(N). Their 
scheme invoked an extra measuring device. As a model device, they chose a precessing 1/2 spin whose precession 
angle should be proportional to the transferred charge. The measurement paradigm is then applied to the device. 



In this way they were able to obtain a satisfactory definition of the full counting statistics Pt(N) with an integer 
number of charges transferred. rThe details of the calculation and a thorough discussion are presented inu 

It was clear to the authors oia that their definition of the FCS does not depend on a specific measurement scheme. 
However, this fact was not explicitly evident. For several years this hindered the impact of these outstanding contri- 
butions. 

One of the authors has recently proposed a slightly different calculation scheme of FCS that does not implicitly 
invoke any measuring device but still leads to the same resultsu. The observation was that the cumulants of the current 
can be obtained as non-linear responses of a system to a fictitious field that can only be defined in the framework 
of the Keldysh diagrammatic technique!! The calculation of FCS can be accomplished with a slight extension of the 
Keldysh technique. This meant some progress since the methods of the Keldysh technique are well elaborated and 
can be readily applied to a variety of physical systems and situations. 

Recent workQ has addressed the charge transfer between two superconductors. The problem can be tackled with 
an extension of the above-mentioned Keldysh technique. The expressions for Pr(N) were obtained. Albeit the 
authors have encountered a significant difficulty with understanding the results in classical terms. The point is that 
the calculation gave negative probabilities. This indicates that the results cannot be interpreted without invoking a 
quantum mechanical description of a detector. 

All this suggests that the quantum mechanical concept of full counting statistics shall be refined and the generality 
of previously used definitions shall be accessed. This is done in the present article. 

To preserve generality, we focus on the full counting statistics of an arbitrary quantum mechanical variable A. Then 
the result does not have to be discrete, and, strictly speaking, no counting takes place. We keep the term "counting" 
for historical reasons. 

We show that introducing a detector is absolutely necessary for defining such statistics. On the other hand, the 
answer does not depend on details of the detector. This allows for a separation of the system measured and the 
measuring device. We develop an exact quantum mechanical description of the system in terms of a path integral 
over detector variables and derive our results from this description. We show that the classical interpretation of full 
counting statistics is only possible in the presence of a certain symmetry. For the full counting statistics of the electric 
current, this symmetry is gauge invariance. In superconductors the gauge invariance is broken and the full counting 
statistics must be interpreted along quantum mechanical lines. 

It is the main message of our paper that this interpretation problem does not make the concept of full counting 
statistics useless and/or unphysical. We show that it is the full counting statistics that completely determines the 
evolution of the density matrix of the detector. We show that thereby the statistics is observable in real experiments. 
We propose and discuss two concrete measuring schemes. 

The paper is organized as follows. We start with a general compact discussion of the interpretation problems. We 



present the model of the detector in section III . We present the approach and derive general relations for the detector 
density matrix in the subsequent section. Section M defines the FCS and presents its interpretation. The subsequent 
sections provide examples of FCS for an equilibrium system in the ground state,a normal conductor and a harmonic 
oscillator. The characterization of the FCS is considered in sections IX and [Xj where two concrete measurement 
schemes are discussed. 



II. GENERAL DISCUSSION 



It is not a priori clear why the operator definition (hi) produces senseless results. We list below possible intuitive 
reasons for this. To start with, the paradigm concerns an instant measurement. The operator definition (yj) is not 
local in time and accumulates information about the quantum state of the system for a (long) interval of time. The 
applicability of the paradigm is therefore not obvious. For instance, the averages of powers of Qt r can be expressed 
in terms of correlators of currents 



< QZ > = / dh...dt N < I(t x )I(t 2 )...I(t N ) > (2) 

Jo 

Usually the causality comes into quantum mechanics via time ordering of operator products. There is no time ordering 
of current operators in (g). This may indicate implicit problems with causality. The second reason is as follows. It 
seems obvious that the time integral of I can be associated with a physical operator of charge. For an arbitrary 
operator A it may be difficult to find such a physical associate. However, the integrals of A can be measured, and 
statistics of those can be accumulated. 

The only way out seems to extend the original quantum mechanical system by adding a "detector". Although 
it is very physical, mathematically this means an eventual doubling of the original Hilbert space by adding extra 



degrees of freedom. In fact, this has been done inH by introducing the spin- 1/2. Below we adopt the simplest way to 
provide such an extensionu. Given a quantum mechanical system, one introduces an extra variable x, whose operator 
x commutes with all operators of the system. We assume that the canonically conjugated variable, q, ([x,q] = i, in 
units with h — 1) can be measured according to the paradigm. Next we introduce an interaction between the system 
and the detector in such a way that in the time interval (0, T) the Heisenberg equation of motion becomes 

hit) = A(t) (3) 

thus simulating Eq. ([[]). In this way we avoid all possible difficulties with misinterpretations of the paradigm. The 
integral of A(t) is now correctly associated with an operator of a value that can be a priori measured. 

Albeit there is a price to pay. As we show below, the FCS can be defined in this way as an operator that relates the 
density matrices of the detector before and after the measurement. In general, it is not the same as the probability 
distribution of the shifts of the detector momentum which may be defined in the classical limit. The FCS can be 
interpreted in such terms only under certain conditions, which are satisfied for the statistics of current in normal 
conductors. However, the definition we adopt is sufficiently constructive since it allows to predict the result of any 
realistic measurement. 

To make the detector more realistic and thus to show the generality of the results, one shall introduce internal 
dynamics of the detector variable. This dynamics would make the detector a non-ideal one: the readings may differ 
from the definitions (nh and (p3). The path integral approach we describe below provides the most convenient way to 
incorporate this internal dynamics. 

III. THE MODEL 

Let us now specify our detection model. The detector will formally be a particle with the Hamiltonian ^— + V(x). 
q and x obey canonical commutation relations. It shall be coupled to the system only for a certain period of time T 
and be decoupled adiabatically for earlier and later times. We introduce a smooth coupling function arit) that takes 
the value 1 in the time interval < t < T and is zero beyond the time interval t\ < t <t%- The values for t\ < and 
£2 > T are chosen in a way that provides an adiabatic switching. The whole Hamiltonian therefore reads 

H(t) = H sys - a T (t)xA +^ + V(x) (4) 

We stress that this detector model is very general. It can easily be applied to describe a variety of electrical circuits 
and mechanical systems. In particular, it may describe selective detectors: by choosing V(x) — ?jtu 2 x 2 we describe 
a harmonic oscillator that selectively detects a signal frequency u). In principle it would also be possible to treat 
detectors that consist of more than one degree of freedom in the formalism we are going to develop. 
The coupling is chosen such that the Heisenberg equation of motion for the detector quasimomentum q 

m = a T (t)A(t)-V'(x(t)) (5) 

reproduces Eq. (||) if the internal dynamics of the detector is disregarded. Eq. (g) suggests that the statistics of 
measurements of the detector's momentum after uncoupling it from the system will tell us about the full counting 
statistics we are interested in. The coupling term can be viewed as a disturbance of the system measured by the 
measuring device. To minimize this disturbance, we would clearly like to concentrate the detector wave function 
around x = 0. The uncertainty principle forbids us, however, to localize it completely. Indeed, in such a way we 
would loose all information about the momentum, which we are going to measure. This is a fundamental limitation 
imposed by quantum mechanics, and we are going to explore its consequences step by step. To concentrate fully on 
the quantum nature of the detector we explore in details the limit of infinite detector mass. In this limit, x = and 
the detector can classically not disturb the system. 

IV. THE APPROACH 

Our goal is to relate the density matrices of the detector before and after the measurement. If there were no system 
to measure we could easily express this relation in the form of a path integral in the (double) variable x{t), like it has 
been done many years ago ircl. The expectation is that the whole information about the system to be measured can 
be compressed into an extra factor in this integral, the so-called influence functional. This functional will depend on 



x(t) rather than on the internal dynamics of the detector. This will make explicit the separation between the detector 
and the measured system. To make contact withu, we would like to present this influence functional as an operator 
expression that involves the system degrees of freedom only. We denote the initial detector density matrix (at t < ti) 
by p ln (x + ,x~) and the final one (at t > £2, after having traced out the system's degrees of freedom) by pf(x + ,x~). 
R is the initial unperturbed density matrix of the system. The entire initial density matrix is assumed to factorize, 
D = Rp m . 

We start out by inserting complete sets of states into the expression for the time development of the density matrix 

pf{x + ,x~) = Tr 

bystcm 

< x+ 1 f cxp{-i \l[ dt [H sys - a T (t)xA +£; + V(x)] } D T~ exp{i j£ dt [H sys - a T (t)xA +£^ + V(x)]}\x~ > (6) 

in the usual manner. ~T(T) denotes (inverse) time ordering. As the complete sets of states we choose product 
states of any complete set of states of the system and alternatingly complete sets of eigenstates of the position or 
the momentum operator of the detector. Those intermediate states allow us to replace the position and momentum 
operators in the time development exponentials by their eigenvalues. We can then do the integrals over the system 
states as well as the momentum integrals and arrive at the expression 

pf(x + ,x-)= f D[x+(t)} D[x-(t)] e^{if 2 dt^[{x + f-{x-f]+V{x + (£))-V{x-{t))} Tr 

J x+(t 2 )=x+ x-(t 2 )=x- Jtx 2 S y stcm 

l^exp{-i J' dt[H sys -a T {t)x+{t)A}} R T~ exp{i [ ' dt[H sys - a T (t)aT (t)A] } p m {x+ (h) 7 x~ {h)) (7) 

We rewrite the expression as 

p f (x + ,x-)= / dx+dx^ K(x+,x-;x+,Xi,T)p m {x+,Xi) (8) 

in terms of the kernel 

K(x+,x-;xf,Xi,T) = J D[x + (t)} D[x~(t)] 

^ x Jr (t2)—x Jr ,x Jr (ti )—x^ x~ (^2)— x~ ,x~ (ti)— a; j 

F[x+(t),x-(t),T] exp{i f 2 dt^-[x+(tf -x-(tf] + V (x + (t)) - V {x~ (t))} (9) 

Jt! 2 

that contains already the desired influence functional 

/•*2 /"*2 

F[x + (t),x-(t),T] = Tr Texp{-i dt[H sys - a T {t)x + (t)A]} RT~ exp{i / dt[H sys - a T (t)x-(t)A) } (10) 

System J t± J tl 

We have therefore accomplished the task: system and detector are presented as separate terms under the sign of path 
integration. 

Seeking further simplifications, we proceed by taking the limit of infinite detector mass and considering a free 
particle, so that V(x) — 0. We find that the kinetic term in Eq. (||) then suppresses all fluctuations in the path 
integral. The details of the calculation needed to establish this are provided in Appendix A. Therefore, the kernel 
K(x + 7 x~ , xf, Xi, T) becomes local in position space, does not depend on the detector and reads 

K(x + ,x-,x+,x^,T) = S(x + -xf) 6(x~ -%) P{x + ,x~,T) (11) 

with 

rt-2 />t2 

P{x+,x-,T)= Tr T exp{-i / dt\H sys - a T {t)x + A] }RT cxp{ i / dt\H sys - a T (t)x~A]} (12) 

Sy stcm Jti Jt! 

If we compare expressions ( pX3| ) and ( ^2| ) we see that the simplification is that the x ± (t) do not depend on time. So the 
detector is somehow "slow" and the variable x(t) quasi-stationary. This gives physical sense to the limit we consider. 
The only requirement is eventually time scale separation between the system and the detector. The typical time scale 



of the internal dynamics of the detector shall exceed that one of the system. Under these conditions, we can apply 
the quasi-stationary approximation and replace F with P. 

It is constructive to rewrite now the density matrices in Wigner representation 



/dz z z 

— e lqz p(x+-,x~-) 



(13) 



and define 



P(x,q,T) = j "|j>* P{x+ Z -,x- Z -,T). (14) 

This gives the following convenient relation 

p f (x,q) = I dq 1 P(x,q-q 1 ,T)p in (x,q l ) (15) 



V. INTERPRETATION OF FCS 

We adopt the relations (|1J), (|4|) and ( |15|) as the definition of the FCS of the variable A. Let us see why. First let us 
suppose that we can treat the detector classically. Then the density matrix of the detector in Wigner representation 
can be interpreted as a classical probability distribution n(a;, q) to be at a certain position x with momentum q. This 

allows for a classical interpretation of P(x, q, T) as a probability to have measured q = J t Q A(t). Indeed, one sees from 
(|l4j) that the final II(x, q) is obtained from the initial one by shifts in q, P(x, q, T) being the probability distribution 
of those shifts. 

Albeit there is a point which hampers such classical interpretation. The density matrix in Wigner representation 
cannot be interpreted as a probability to have a certain position and momentum or any other probability, since it is not 
positive. Concrete calculations given below illustrate that P(x, q, T) does not have to be positive as well. Consequently, 
it cannot be interpreted as a probability distribution. We stress that the absence of a classical interpretation does 
not hamper a "predictive power" of FCS: having it in hands, one can predict the results of measurements using Eq. 
(O). The reverse problem is to extract the FCS from the results of measurements. We discuss two schemes for this 
in sections IX and |x|. 



There is an important case where FCS can indeed be interpreted as a probability distribution. In this case P(x, q, T) 
does not depend on x, P(x,q,T) = P(q,T). Then, integrating Eq. (119) over x, we find 



Uf(q) = Jdq'P(q-q',T)Il in (q') (16) 

with 11(g) = J dxp(x, q). Therefore, the FCS is in this special case the kernel that relates the probability distributions 
of the detector momentum before and after the measurement, H m (q) and H^(q), to each other. Those distributions 
are positive and so will the FCS P(q,T). Note, that a generic case in which this interpretation is appropriate is the 
measurement of a static operator A, that is, [H ays , A] = 0. Then the FCS is obviously x-independent and it can be 
interpreted as the probability distribution of outcomes of an instant measurement of A. 

Since we study FCS of a stationary system and the measurement time T exceeds timescales associated with the 
system, the operator expression in Eq. (O) can be seen as a product of terms corresponding to time intervals. 
Therefore in this limit of T — > oo the dependence on the measuring time can be reconciled into 

P(x + ,x-) = e- e ( x+ ' x ~ )T (17) 

where the expression in the exponent is supposed to be big. Then the integral ( |l4|) that defines FCS can be done in 
the saddle point approximation. If one defines A — q/T, that is, A is the result of measurement of J A(t)dt/T, the 
FCS can be recast into the similar from 

P(x,A,T) = e- iix ^ T (18) 

where £ is defined as a (complex) extremum with respect to (complex) z: 

£ = e^tr{£(x+^,x- |) - iAz} (19) 

z 2 2 



The average value of A and its variance (noise) can be expressed in terms of derivatives of £ : 

,. d£(x + z/2,x-z/2) a2 d 2 £{x + z/2,x~z/2) 

<A>=hm i '-2 L^ ; T<CA 2 >=hm 1 ^ '—!-. (20) 

2^0 idz 2^0 dz 2 

Indeed, the quantity P{x + , x~ , T) is the generating function of moments of q. It is interesting to note that in general 
this function may generate a variety of moments that differ in the time order of operators involved, for instance, 

xM xN-M pT 

Q N m = {~i) N Jm d{x _ )M d{x+)N _ M P{x + ,x-) = J dh...dt N < T-{A{t 1 )...A{t M )}-f{A{t M+1 )...A{t N )} > (21) 

The moments of q are expressed through these moments and binomial coefficients 

iNrtN 



Q W = f dqq N P(0, q, T) = 2- N Y, CM 1 22 i 

J M 



VI. FCS OF AN EQUILIBRIUM SYSTEM IN THE GROUND STATE 

To acquire a better understanding of the general relations obtained we consider now an important special case. We 
will assume that the system considered is in its ground state \g >, so that its initial density matrix is R = \g >< g\. 
This assumption allows for an easy calculation of the FCS. We have supposed that the coupling between the system 
and the detector arises adiabatically. Then the time development operators in ( |l2| ) during the time interval 1 1 < t < 
adiabatically transfer the system from \g > into the ground state \g(x ) > of the new Hamiltonian -ff sys — x A . In 
the time interval < t < T the time evolution of the resulting state has thus the simple form 

exp{-i t {H sys + x ± A)} 1.9^) > = e~ UE{x± ^ ^(x*) > . (23) 

E(x ) are the energies corresponding to \g(x ) >. This provides the main contribution to the FCS if the measurement 
time is large and the phase acquired during the switching of the interaction can be neglected in comparison with this 
contribution. Therefore we obtain 

P(x+,x-,T) = e -iT(B(. + )-E(0). (24) 



First we observe that P(x,q,T) is a real function in this case, since the exponent in (25) is anti-symmetric in z. A 
first requirement for being able to interpret P(0, q, T) as a probability distribution is therefore fulfilled. However, 
the same asymmetry assures that all even cumulative moments of q are identically zero, whereas the odd ones need 
not. On one hand, since the second moment corresponds to the noise and the ground state cannot provide any, this 
makes sense. On the other hand, this situation would be impossible if P(0,q,T) were a true, positive probability 
distribution. _ 

Belzig and NazarovQ encountered this situation analysing the FCS of a general super-conducting junction. In a 
certain limit the junction becomes a simple Josephson junction in its ground state. In this mere limit the interpretation 
of FCS as a probability distribution does not work any longer. Fortunately enough, the relation ([15|) allows us to 
interpret the results obtained. 

To keep simple things simple, we assume the function E(x) to be analytic and expand it in its Taylor series. We 
also rescale q as above, q — A/T. For the FCS we therefore have 

P(X,A,T)= f^l e "AT . e - l T(E'( x ) z +E'" (x )^/2A+...) (25) 

If we take the limit T — ► oo now, we note that all the terms involving higher derivatives of E(x) are negligible and 
arrive at 

lim P(x,A,T) =5{A-E'{x)) (26) 

T — >oo 

According to the Hellman-Feynman theorem E'(x) =< g(x)\A\g(x) >. As one would expect, in this limit the 
measurement gives the expectation value of the operator A in a ground state of the system that is somewhat altered 
by its interaction with the detector at position x. Therefore the resulting dispersion of A will be determined by the 
initial quantum mechanical spread of the detector wave function. Consequently, the error of the measurement comes 
from the interaction with the detector rather than from the intrinsic noise of the system measured. 



VII. FCS OF THE ELECTRIC CURRENT IN A NORMAL CONDUCTOR 

A complementary example is a normal conductor biased at finite voltage. This is a stationary non- equilibrium 
system far from being in its ground state. Here we do not intend-to go to microscopic details of the derivation. Our 
immediate aim is to make contact with the approaches of Refs.Qti. We keep the original notations of the references 
wherever it is possible. 

The starting points of the approaches differ much. Levitov and Lesovik propose a detector model where the z- 
component of a spin-1/2 creates a local vector potential which is felt by the electrons. This corresponds to a total 
Hamiltonian of the form 

H = H sys — —<r z I 

which is studied at different coupling constants A. The referenced starts with an extension of the Kcldysh technique 
to formally unphysical systems where the evolution of the wave function in different time directions is governed by 
two different Hamiltonians 

#± = H sys ± X I (27) 

and shows that the so defined Green functions can be used to generate moments of /. This is to be compared with 
our detection model. 

Despite different starting points, all three approaches quickly concentrate on the calculation of the quantity 

< exp(-zFiT) exp{iH 2 T) > (28) 

which is identical to our definition of P(x + , x~) provided that iJ 12 = H sys — x^I. This quantity is denoted by xW 
iria and by exp(— S(x)) ina- It is used to define the FCS and we see now that the final result does not depend on the 
starting point. n 

To give a concrete expression for FCS, we just make use of Eq. (37) oD. We consider the FCS of the current in 
a multi-mode constriction which is characterized by a set of transmission coefficients T n . In general, the answer is 
expressed in terms of energy-dependent electron filling factors Uril) on the right (left) side of the constriction 

P(x+,x-,T) = exp{— V / dE\n(l + T n (e^-* + y e -l)n R (l-n L ) + T n (e^ + - x ^/ e -l)n L (l-n R ))} (29) 

We stress that this expression depends on x + — x~ only. This is a direct consequence of gauge invariance. Indeed, 
in each of the Hamiltonians the coupling term is the coupling to a vector potential localized in a certain cross-section 
of the constriction. The gauge transform that shifts the phase of the wave functions by x if on the right side of 
the constriction, eliminates this coupling term. This transform was explicitly implemented intl. Since there are two 
Hamiltonians in the expression, the coupling terms cannot be eliminated simultaneously provided x + =/= x~ . However, 
the gauge transform with the phase shift (x + + x~)/2 makes the coupling terms depending on x + — x~ only. 

Since P(x + ,x~,T) depends on x + — x~ only, the FCS P(x,q,T) does not depend on x. As we have stressed in 
section M, this enables one to interpret the FCS as a probability distribution. 

Superconductivity breaks gauge invariance, thus making such an interpretation impossible. 

VIII. FCS OF A HARMONIC OSCILLATOR 

Let us now illustrate the measuring process proposed with a simple example. We consider the measurement of the 
position of a harmonic oscillator in its ground state. That is, Hq — ||y + ^Mlo 2 X 2 is the Hamiltonian of the system 

and A = X shall be measured. The measurement is effected by coupling the oscillator to a detector with position 
variable x during a time interval of length T. In the limit of infinite detector mass the entire Hamiltonian in the 
measurement period reads then 

Q 2 1 

H = — + -Muj 2 X 2 + xX. (30) 

2M 2 v ; 



The perturbed ground state \g(x) > is in this simple example obtained by shifting the original ground state wave 

l 

2ffi? 



function by x/Mcu 2 in A-representation. Its energy is E g {x) = E g (0) — 2J J . x 2 . We then find from (p5|), that 



P{x,q,T) = 8{q + xT/ML0 2 ) (31) 

Following our first classical interpretation of P(0, q, T) we would now conclude, that a harmonic oscillator in its ground 
state does not transmit any fluctuations of its position variable to the detector and that the detector's wave function 
is not altered by the oscillator. Let us, however, calculate the read off of the detector after the measurement just 
described, i.e., its final momentum distribution. We have to integrate the final density matrix found from (n5|) over 
x. As the initial detector state we choose a Gaussian wave with uncertainty Ag in the momentum: 

p in (x, q) = exp[-^2 - 2(A<?) V] (32) 

The final momentum distribution 

n/( g ) = exp[ 2(Ag)2 + r - g 2M2 ^ 4(Ag)2 ] (33) 

is again a Gaussian. Its width Aq* , however, increases in time, for big times linearly with T: 

^) 2 = ^ 2 + 4mJ W (34) 

This relation clearly demonstrates the relevance of the detector's influence on the system due to its quantum nature. 
The initial detector momentum is known only with an uncertainty Aq. This uncertainty is a source of measurement 
error and is described by the first term of Eq. (Q). We would clearly like to minimize this error by choosing Aq 
to be small. Doing so, we increase, however, the uncertainty in the detector position Ax > l/2Aq. According to 
Eq. (g^), such a spread in the detector position will disturb the system. Since the oscillator is in its ground state 
this disturbance will drive it into ground states of new Hamiltonians H$ + xX. For every detector influence x we 
will therefore measure a different expectation value E'(x) —< g{x)\X\g{x) >. The read off of the detector will be a 
superposition of measurement outcomes corresponding to all those different influences. The resulting broadening of 
the final detector momentum distribution (Aq^)^ x w Ax E"(0)T is accounted for by the second term in (|34|). We 
conclude that the the quantum fluctuations of the detector set an upper bound on the accuracy of the measurement 
process. The actual value of this bound depends of course on the system measured. It will vanish if the FCS is 
x-independent and a classical interpretation of the process is appropriate. 

We have seen that the FCS is at least not directly observable in a single experiment. It is buried under some 
additional broadening which we could fully understand in terms of the disturbance of the system by the detector. It 
therefore does not tell us anything about fluctuations of the observable A in our system, the quantity we want to 
measure after all. In the next section we will therefore try to get rid of this broadening and measure the function 
P(x,q,T) itself, which is more specific to the system. A certain deconvolution procedure will be developed. 

IX. CHARACTERIZATION OF THE FCS; FIRST SCHEME 

As we have already seen, the full counting statistics P(x, q, T) proposed above allows us to predict the outcomes 
of measurements with a quite general detector and resolves inconsistencies that arose in earlier interpretations. It 
remains to show now, that it is a physical, real object in the sense that it is experimentally observable. 

We would like to suggest two schemes for measuring the FCS. We start from relation ( |15J ) between the initial and 
the final density matrix. Writing this equation in (x + , ai~)-space, we find the simple expression 

P(x+,x-,T)= P f n { * + + > X l\ (35) 

p m (x^,x ) 

or 

P(^T)=[pe^ if + Z ^ X - z/ % (36) 

J 2tt p m (x + z/2,x — z/2) 

for the kernel. We would already be done if we could measure all elements of the detector's final and initial density 
matrices. This is not possible in general, however. By successively measuring a certain observable we can measure the 
diagonal elements of the density matrix in a basis of eigenstates of that observable, but not the off-diagonal entries. 
We can therefore measure the functions II(<7), but not p itself. 



The key idea that we will pursue to solve this problem is to repeat the same measurement many times for shifted but 
otherwise identical initial detector density matrices. We suggest to repeat the measurement of the final momentum 
distribution II ^ (q) for a number of initial density matrices that differ only in the expectation value xq of the position 
of the detecting particle and call it 

Tf(x , q, T) = n£(c?, x ) = Jdx dq' P(x, q - q' , T) p m (x ~ x ,q'). (37) 

This way we expose the system during the measurement to different detector influences and one can hope that by 
doing so this influence can be identified and eliminated by some kind of a deconvolution procedure. Defining the 
Fourier transform of T' (xo, q, T) with respect to both of its variables 

f '(go, z,T) = ±-f dx Q dq e ix °«°- iz « r' (x , q, T) (38) 

we find, that the FCS can indeed be reconstructed from this function by means of the relation 



P(x, q,T) = ±- [ dq dz e*"-**" T ™' Z '\ ) (39) 

2tt J p m (q ,z) 



where 



p(q Q ,z)^jdxe^ x p(x + ^x-^). (40) 

To interpret the result of the measurement, we still have to know the full initial density matrix of the detector. This 
should be feasible, however. One might cither prepare the detector initially in a specific, well-known state, or one might 
let the detector equilibrate with an environment. Then the initial density matrix is stationary, = [p m , H] oc [p m , q 2 }, 
it will be diagonal in a basis of momentum eigenstates and can be determined by a momentum measurement only. 
This proves that the FCS really is an observable. 

To illustrate the procedure we want to apply it now to our example of the harmonic oscillator already introduced 
earlier. Like there IP (q) we calculate now the function F* (xo, <?, T) that would be measured in the actual experiment. 
We shift the initial oscillator wavefunction by x$ in x-space and find that 

rff t\ r -(g + xpT/Mio 2 ) 2 1 

r/(, , g ,r)^exp[ (2(Ag)2 + r2/2M2 ^ 4(Ag)2) ]. (4i) 

On transforming this into Fourier space it becomes 

ff(q ,z,T) = exp[-^^ - 8 M^(A g )^ 6 ^+Tz/Muj 2 ). (42) 

Employing now Eq. ( p9[) with p m (q ,z) = exp[— ' 2 slAg)^ ] we m£ ieed recover the desired FCS Eq. (J31|). 

X. SECOND SCHEME 

We see also a possibility of measuring P(x, q, T) without knowledge of the initial state of the detector. We can 
apply it, however, only to the special case in which the system is in one state only, for example its ground state, or in 
a mixture of a limited number of discrete states. This condition may be satisfied at low temperatures or if we have 
dissipation in the system which drives it into the ground state. We explain the idea for the case of a system in its 
ground state. Then we have the explicit expression ( |24| ) for the time evolution and we find, that 

r-f (xo, q, T)= I dx dz dq' e -i<1-l')-iT{E'{x) z+±E'"(x) z*+...) p in^, _ ^ ^ (43) 



r*(xo, q,T) again being the final momentum distribution for shifted initial detector wave functions. In the limit of 
big T we find with @ 

lim T f (x , q, T) oc / dx p m (x -x ,q- TE'(x)). (44) 



This formula suggests that we can measure the function E'(x) arbitrarily exact in the limit of long measurement time 
T by determining the peak of the final momentum distribution. The only assumption we have to make about the 
initial detector density matrix now is, that it is well centered around x = and that it falls off sufficiently fast for 
momenta higher than some arbitrary Aq. We want p m to be peaked in x-space such that E'(x) is measured at the 
point xq only (we assume E(x) to be analytic). Of course, this means, that the width Aq in momentum space of p m 
and therefore also of p T will be rather wide. This is not a problem, however, since the position of the peak that we 
can measure with uncertainty Aq increases linearly in time. The uncertainty in our knowledge of E'(x) is therefore 
determined by-j? and tends to zero in the infinite time limit. 

Integrating E'(x) we can then reconstruct the FCS for arbitrary detection times 



f dz i-x+z/2 

P(x,q,T)= / — exp{iqz-iT / dx' E'(x')} 

J 2n Jx-z/2 



s/2 

If the system is in a mixture of N distinct states, the expression for P(x, q, T) will be a sum of terms of the form 
( J24| ) with different functions E n (x). We will in general find N distinct peaks in the final momentum distribution 
allowing us to measure all N functions E' n (x). Again, we can reconstruct P(x,q,T) for arbitrary T. 

XL CONCLUSIONS 

In this article we have been studying the process of measuring a general quantum mechanical variable A by means 
of a detector that is coupled to it for a certain period of time. We have been focusing on a simple model describing a 
detector without internal dynamics which is nevertheless applicable to a variety of experimental situations. Moreover, 
the formalism that we have presented is general enough to allow for the description of a generic detector. We have 
predicted the outcome of such a measurement. This prediction involves a new object that we want to call full counting 
statistics (FCS) of the variable A and that is an extension of another function proposed earlier in this context. This 
extension basically consists of accounting for the detector influence on the measured system. We find that the interplay 
of this influence with the quantum nature of the detector hampers in general a classical interpretation of the detector 
read-off. This way we have been able to remove inconsistencies, to be more specific, negative probabilities, that arose 
in the previous interpretation. Finally, we have shown, that this FCS is not only a theoretical construct that helps to 
predict results of certain measurements, but that it is an observable itself. It can be measured in experiments. 
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APPENDIX A: 

In this appendix we give a detailed derivation of the infinite mass limit of the kernel (J9|). 
First, we define a Fourier transformed influence functional 

F[k + (jt) t k-{t),T\= J D[x+(t)]D[x-(t)}F[x + (t),x-(t),T] 



exp{i / dt[x + {t)k + (t)~x-(t)k-(t)]} (Al) 

and correspondingly 
F(q + , q~ , Qi, qi , T) = / dx + dx~ dxj dx^ exp{ix + q + — ix~ q~ — ix^q^[ + iXiq^}K(x + , x~ , x[ , xj, T) (A2) 
The k (t) are general functions on the interval [ti,^]- Inserting the identity 
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exp{z J' 2 dt | x{tf) = J D[q(t)) exp{* ^ dt [-^- - q(t)i(t)} } (A3) 

and using (0) we derive then 

K(q+,q-,q+,q-,T)= f D[q+(t)} D[q~(t)] 

J q +(t 2 )=q+, q +(t 1 )=q+ q-(t 2 )=q-.,q-{t 1 )=q 1 

F[q+(t),q-(t),T] expjz f t2 dt(&--^-)} (A4) 

J tl 2m 2m 

In the infinite mass limit, the kinetic term in this expression will disappear. Now, we change integration variables 
from -D^*] to -D^] and call fc ± = q . Then, 

K(q+,q-,q+,q^,T)= f D[k+ (t)] D[k~ "(<)] F[k+(t), k~(t), T] (A5) 

"" p 2 dtk+{t)=q+-q+ f* 2 dtk-{t)=q--q- 

We can represent the functions fc ± (t) by their Fourier series, fc ± (t) = z~2^Lo k^ cos n * t _* - Changing the inte- 
gration variables in (A5) to the coefficients k^ in this expansion, we notice, that only the integrals over the zeroth 



components k are constrained by the boundary conditions. We can therefore do all the integrals over higher Fourier 
modes in (|9|) and obtain 

K{q + ,q-,qt,qi,T)= f D[k+] D[fc~] D[x+] D[x~] 
F[x+ft),^(t),r]exp{ t(t2 ~ tl) f;(fc+x+-fc- a: -)+,[4( g +- g +)-x -( g -- gl -)]} (A6) 



n=l 






We have also expanded the functions x (t) in Fourier series now, x (t) = z~2^Lo x n cos 

We see, that the /^-integrations result in ^-functions that constrain the x^, n ^ 0, to zero and allow us to do the 
corresponding x^ - integrals: 

K{q + 1 q~,qi,q^ 1 T)= dx£ dx cxp{i[a;+(q + - qf ) - Xq (q~ - g7/ )] } 

Tr T cxp{-i / dt\H sys - a T (t)x£A] } R T exp{i / <ft[ff sys - a T (t)x^ A] } (A7) 

s y stom Jtt Jtt 

or 

if(x + ,a;-,a;^,a;r,T) = 6(£ + - a^ ) <J(aT - »7/) P(x + ,x~,T) (A8) 



with 



/■*2 _ /"i2 

P(a: + ,ar,T) = Tr T exp{-i dt[H sys - a T {t)x + A]} RT exp{ i/ dt[F sys - a T (t)x _ i]}. 
system Jt! Jti 

Thus, the claimed locality of the kernel K (x + , x~ , x^ , x^ ,T) is established. 



(A9) 
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